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Abstract 

The inflationary cosmology is analyzed from the point of view of squeezed 
quantum states. As noted by Grishchuk and Sidorov, the amplification of 
quantum fluctuations into macroscopic perturbations which occurs during 
cosmic inflation is a process of quantum squeezing. We carefully develop the 
squeezed state formalism and derive the equations that govern the evolution 
of a gaussian initial state. We derive the power spectrum of density pertur- 
bations for a simple inflationary model and discuss its features. We conclude 
that the squeezed state formalism provides an interesting framework within 
which to study the amplification process, but, in disagreement with the claims 
of Grishchuk and Sidorov , that it does not provide us with any new physical 
results. 
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1 Introduction 



One of the impressive features of an inflationary cosmology is the prediction 
of a set of perturbations on the background Robert son- Walker metric. These 
perturbations are produced via the amplification of ground state quantum 
fluctuations during the inflationary period. This process has been widely 
studied and there is broad agreement regarding both methods and results 
[|]. The actual perturbations predicted depend on details of the inflationary 
period. A cosmology with a period of simple exponential inflation and with 
cold dark matter forms the basis of the "standard CDM" model for the 
formation of galaxies and other structure in the universe. This model has 
enjoyed great popularity, but it is also coming under increasing pressure from 
astronomical observations fU , II , [II , 11 ■ 

Recent work by Grishchuk and Sidorov || , [0] has suggested that impor- 
tant quantum effects have been neglected in the standard approach. These 
authors claim that, because of quantum squeezing, inflation predicts features 
in the perturbations which have not been properly taken into account and 
which could result in striking observational consequences. In particular, they 
emphasize the phenomenon of desqueezing which leads to approximate zeros 
in the power spectrum at calculable wavelengths. 

We have systematically investigated the inflationary cosmology from the 
point of view of quantum squeezing, using Bardeen's gauge invariant variables 
[§]. We have found that indeed each mode of the perturbed field evolves as a 
squeezed state during the inflationary period but that the features discussed 
by Grishchuk and Sidorov in |3j and J7[ are well known ones, which are 
essentially classical in nature. Although we note in section 6 an isolated error 
in the literature which may have prompted much of Grishchuk's criticism, we 
argue that the error can be (and usually is) avoided without appealing to the 
formalism of squeezed quantum states. We conclude that this perspective 
offers nothing more then an alternative set of words (and variables) with 
which to discuss the inflationary universe. We do however find the squeezed 
state formalism well suited to the problem || and it may prove useful in 
future work. 

The structure of the paper is as follows. In section 2 we look at a sim- 
ple mechanical system - the inverted harmonic oscillator - and show how 
it exhibits squeezing behavior at both the classical and quantum levels. In 
sections 3 and 4 we use the formalism of gauge invariant cosmological pertur- 
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bations as presented in [|I(J , to construct the Hamiltonian operator. We then 
set up the time evolution operator and show that it can be factorized into 
a product of a squeeze operator and a rotation operator which are charac- 
terized in terms of the squeeze factor Rt, squeeze phase and the rotation 
angle 6g. R^ gives us a measure of the excitation of the state while $g gives 
us a measure of how the excitation is shared between canonical variables. We 
show how the evolution of the state can be characterized by a set of coupled 
first order ordinary differential equations for and 0£. In section 5 

we study the behavior of this system of ODE's, identifying different regimes 
according to the scale of the perturbations: on scales larger then the Hubble 
radius the squeeze phase freezes out and the squeeze factor grows; on scales 
smaller than the Hubble radius the squeeze parameters oscillate. 

Having gained some insight into what to expect generically in such mod- 
els we look at a simple inflationary model with baryonic matter coupled to 
photons (without dark matter) such that the evolution of perturbations can 
be well approximated by a single collective scalar field. We generate some 
typical power spectra, |5fc| 2 , and see that they are Harrison-Zeldovich on su- 
perhorizon scales (|<5fc| 2 cx k - no oscillations) and exhibit standard sound 
wave oscillations on subhorizon scales. 

In section 6 we discuss the desqueezing effect emphasized by Grishchuk 
and Sidorov and argue that it a familiar one properly taken into account in 
standard calculations. In section 7 we attempt to clarify the claim that these 
effects are of a distinctly quantum mechanical origin. We comment, using the 
language of squeezed states, on the classicality of the harmonic oscillator; we 
note that, for large squeezing, the squeezed state satisfies the WKB criterion 
for classicality. This is equivalent to the WKB classicality at late times in an 



inverted harmonic oscillator studied by Guth and Pi in Ref. |TTJ. The point 
of this section is to explain that the apparently very quantum mechanical 
squeezed state is in fact classical in the sense with which cosmologists are 
familiar. That the truly quantum mechanical features of these states which 
are probed, for example, in quantum optics might have cosmological impli- 
cations is a fascinating claim but one which has no substance at present. In 
section 8 we summarize briefly and conclude. 
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2 The Single Inverted Harmonic Oscillator 



The aim of this section is to familiarize the reader with the language of 
squeezed states. We apply the squeezed state formalism to a simple system - 
the inverted harmonic oscillator. We will show first how this system exhibits 
squeezing behavior at the classical level. We show how this behavior is due 
to the presence of one growing and one decaying solution and that essentially 
the same behavior carries over to the quantum mechanical system. 



2.1 Classical 

The inverted harmonic oscillator (with unit mass and spring constant) is 
described by the Hamiltonian: 

H= P -- q -. (1) 

2 2 K J 

A convenient choice of variables is 

b ± = -j=(p±q). (2) 

The general solutions are: 

b + (t) = M0)e* &-(*) = 6-(0)e"*. (3) 

The evolution of the inverted harmonic oscillator is illustrated in Fig 1, which 
shows the trajectories in phase space of a few representative solutions. The 
phase space can be labeled equally well by p and q or b + and 6_ (the rotated 
axes). As time goes on the value of b + gets exponentially large, while the 
value of 6_ gets exponentially small. This is because all (but one) of the 
solutions eventually go "rolling down the hill" . As this occurs, p and q each 
grow exponentially, while their difference exponentially approaches zero. 

The trajectories in Fig 1 describe squeezing in the sense that they get 
closer together in the b- direction and further apart in the b + direction. For 
example, the circle in Fig 1 evolves into the squeezed shape above it after a 
period of time. Any probability distribution in phase space will eventually 
become squeezed along the p — q axis as the system evolves. 
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2.2 Quantum 

Now consider the quantum system described by Eq. [I]. Using the usual a 
and a) defined for the right-side-up harmonic oscillator we find that 

H = P -- q - (4) 

= ^(av^-k.), (5) 

We have written the Hamiltonian in this way because this is the form di- 
rectly comparable with the more general squeeze Hamiltonian which we will 
consider. 

If the system starts in the vacuum state annihilated by a (which is just 
the gaussian ground state of the right-side-up oscillator) it evolves into a 
"squeezed state" given by 

|tf(t)) = S\0) = e^ 2e ' 2 ^- h ^\0). (6) 

The "squeeze operator" S is specified by two parameters: r, the "squeeze 
factor" , and 0, the "squeeze phase" . For a general squeeze operator r and 
can be complicated functions of time, but in this simple case they reduce to 
r = t and = — 7r/4. 

We now discuss the squeezed state in connection with the Heisenberg 
uncertainty relationship. Using the relation 

<S' aS = a cosh r — a)e 2l< ^ sinh r, (7) 

it can be easily shown that 

p|*(t)> = a(r,0)gfoKt)>, (8) 



where 



It then follows that 



. cosh r + e 2 ^ sinh r 
cosh r — e 2 " 1 ^ sinh r 



[^(t)\p^(t)) = \a(r,<P)\ 2 ^(t)\q 2 \^(t)) (10) 
= — (cosh 2r + sinh 2r cos 20), (11) 
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and the uncertainty relationship is 



({V(t)\q 2 \V(t)}{y(t)\p 2 \m}) 12 = |(l + sin 2 20sinh 2 2r)i (12) 

Thus AqAp ~ ^he 2t for t ^> 1. The initial minimum uncertainty gaussian 
state which "sits at the top of the hill" spreads rapidly in q and p. 
Consider however 

(^(t)|(pcos0-gsin0) 2 |^(t)) = ^e 2r (13) 

and 

(^(t)|(psin0 + gcos0) 2 |^(t)) = ^e~ 2r . (14) 

For <p = — 7r/4 these are just (Ab + ) 2 and (A6_) 2 . Thus in the p — q plane we 
say that the state is squeezed along an axis with slope tan (f>. The fluctuations 
normal to this axis are exponentially small. This behavior mirrors that of 
phase space trajectories for the classical system (see Fig. 1) and likewise 
corresponds to the existence of one decaying and one growing solution. 

The state can in fact be represented as a phase space density, using the 
Wigner function |12j| , for which the contours are ellipses with one axis of 
length e 2r defined by the angle and the other axis of length e~ 2r as in Fig. 
1. The squeezed states which we will consider will have a time dependent <p 
so they can be pictured as ellipses rotating in the phase space. 

Quantum squeezed states generate considerable interest in various areas 
of physics, e.g. nonlinear optics |TJ|], gravity waves |TJ[], flIE| , gravity 



wave detectors, and quantum cosmology [pL7|| . Their striking feature is that 
they exhibit dramatically the Heisenberg uncertainty relation, by allowing 
one variable to have arbitrarily small uncertainty. The conjugate variable 
has a compensating large uncertainty so the Heisenberg uncertainty relation 
is obeyed as an equality. In this sense squeezed states are very quantum 
mechanical. We will discuss the issue of classical vs. quantum aspects in 
more detail in section 7. 
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3 Formalism For Cosmological Perturbations 



The gauge invariant formalism of cosmological perturbations is well suited 
to the study of the evolution of vacuum fluctuations. As discussed in [|K] 



the problem is reduced to the analysis of the evolution of a scalar field with 
a time dependent mass. 

If one looks solely at the scalar degrees of freedom of the metric pertur- 
bations 

^-^)(% 2 -(^._ £|ij) ). (15) 

it is possible to combine the functions <p, ip, E, B into two gauge invariant 
quantities (invariant under local coordinate transformations) 

= -0 — H{B — E') , $ = (f>+(l/a)[(B - E')a}' , (16) 

where 7i = a' /a is the conformal Hubble parameter, a denotes the scale 
factor and ' = d/dr] denotes the derivative with respect to conformal time. 
We can do the same thing with the matter fields; for example with a scalar 
field, (p(x,r)) = (po(rj) + 5<p(x,r)), we can build a gauge invariant quantity 

6<pW = 5(p + <p' (B - E') . (17) 

These gauge invariant quantities can be combined into a single scalar field 

u = aOM2t + z*), (18) 

where 5(p^tt denotes a generic matter field perturbation, z is given by 

z = (a/c s H)[\{H 2 -H')] 1 l\ (19) 

and c s = (Spo/5eo) 1 ^ 2 denotes the speed of sound (in inflation the correct 
equations are obtained by setting c s = 1). The action for the perturbations 
can then be written as 

Spert = \J d A x[{v'f - C 2 s {v,f + (20) 
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which is the action for a free scalar field v with a time dependent mass 
(m 2 = —z"/z) |10| . Up to a total derivative term this action is equivalent to 
the action 

S' pcTt = \ J d'x[{v'f - cl{v,f - 2 Z - z w' + (^)V] , (21) 

which we will find more convenient to work with. We can now proceed with 
the standard quantization. Constructing the Hamiltonian we get 

H =^J d 3 x[n 2 + c 2 s {v,i) 2 + 2^-v7i]. (22) 

Promoting the fields to operators and taking the Fourier decomposition so 
that 

:v t e lk x 



(2tt) 3 /2 "* l 

* - <*») 

we get the two-mode Hamiltonian 

z' 

n k = ^-k^k + c l k2 v- A + ~ + ■ ( 24 ) 

We want to work in the Schrodinger picture, in which the operators v? 
and 7r^ are fixed at an initial time. We define modes with initial frequency 
equal to k which, suitably normalized, give 

% = -»Oo S -oi s ). (25) 



The two-mode Hamiltonian operator can be written in the simple form 
Kg = Hf + Uf = n^ajas + a^a_ n + 1) + <A s ( e - ai ^a s o_ jE - h.c.) , (26) 
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where 




(27) 



Eqs. (BB) and ( B7[ ) describe the generic momentum conserving quadratic 
Hamiltonian for a scalar field. It has a free evolution piece, 7~t^ with a 

time dependent frequency Qj:, and a squeezing piece, TCr , with a coupling 
strength A^(t). The evolution operator produced by this Hamiltonian can be 
factorized in the following way 



where TZ is the two-mode rotation operator defined as 

K[G%] = exp[-i0g(ata^ + a^a_g + 1)] 

and S is the two-mode squeeze operator defined as 

R 



f (e" 2 ^a_^ - h.c. 



(28) 



(29) 



(30) 



This simple decomposition of the evolution operator is a general property of 
momentum preserving quadratic Hamiltonians fL8|. The rotation operator 



alone gives ordinary oscillations (points in the phase space of a classical 
harmonic oscillator rotate about the origin). The squeeze operator alone 
produces squeezing as discussed in Sect. |]. The complete solution to the 
problem we are considering reduces to finding and 9^ as functions of 

time. (Note that is not the Bardeen variable which we shall write $ B !) 



4 Evolution Equations 

In this section we address the generation of cosmological perturbations by 
studying the evolution of the initial vacuum state with the Hamiltonian dis- 
cussed in the previous section. 
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4.1 The Squeezed Vacuum State 

To begin with we have to define the initial conditions of our quantum field 
theory. We assume that all the modes of interest (i.e. the modes on subhori- 
zon scales today) are well within the horizon at the initial time. In this case 
we have k\rj\ 3> 1, which (with c s = 1) implies Qr ~ k ^> Ag ~ (1/|7/|), and 
Eq. (|26|) reduces to the free Hamiltonian ' ■ We then choose for the initial 
state the ground state of the free Hamiltonian, i.e. the Poincare invariant 
vacuum state, which is defined by 

ag|0>in = 0, V*. 

The action of the rotation operator TZ produces an irrelevant phase 

ft[e£]|0> in = e ie S|0) in , 

but, when acted upon by S[R^, the vacuum state transforms into a two- 
mode squeezed state [T8| 

OO 1 

\SS U ) =S[R u ,* u ]\oy m = E— r^(-e 2i ^t^hR n r\n,k;n,-k) , (31) 
where 

OO 1 

\n,k ]n ,-k) = J2-M a -k) n \°^ (32) 

n=0 n - 

is the two-mode occupation number state. This part of the evolution oper- 
ator is responsible for the amplification of the initial vacuum fluctuations; 
momentum conserving pairs of quanta are created. The squeeze factor is 
related to the mean number of quanta, nr, in the squeezed vacuum state 
through the relation 

nj;=(SS\Nj;\SS) = S inh 2 R n . 

4.2 Evolution Equations 

The problem is to determine the functions R^, $g and 0r. The time evolution 
operator is given by the time ordered exponential 

U( V , Vo ) = T expH f dv"H n (r)')) 
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Texp 



exp 



>lo 



a k a -k ~ n - C ' 



I ID 



i I dr]'n % {r]'){aia k 



a -k a -k 



(33) 



We divide the evolution into infinitesimal time intervals [] e. The composite 
property of the evolution operator implies 



U{j] + e, r] ) =U(r) + e, r\)U (77, 7? ) 



(34) 



We can recast this in terms of the squeeze operator S and the rotation 
operator 1Z in the following form 



5[%$ s ]w[e s 

Taking account of Eq. 



(35) 



we infer that for small e: 5R^ ~ 



VLj:{r])e and 5(pj: ~ Using the computation properties of the squeeze and 
rotation operators, the right hand side of (|35| ) can be written as 

RHS = 5(5%, <%] S[Rl 4>\ ~ %] K[el + <%] . (36) 

In order to express the product of the two squeeze operators in terms of a 
single squeeze operator we use the standard composition property, as given 
in 



18 



S[SRg, %] S[R% <j>\ - 56^} = <S[% * s ] n[9 



k\ 1 



where 



e 1 ^ cosh 



cosh cosh SR^ 



+e vv * 



?s) sinh i#sinh tSiZr 



sinh Rt 



sinh i?S cosh 5Rj; 



^ sinh 5ik cosh (37) 



For sufficiently small e we can expand the LHS in SR^ and 58^ to obtain the 
recursion relations 



©jffa + c) 



+ A^Je cos 2(^(77) - $^)) 

- fi ^) e + ^fa) k^Rdv) + coth^(r/)) sin2(^) - <%(r/)) 



1 Note that in contrast to the evolution operator which describes parametric amplifi- 
cation in p8[ , the time ordering problem is non-trivial. 



(3* 
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where Gg = 6^ + 59 'g + 6. The differential form of these equations is 



These are the equations of motion of our system. The analogous equations 
for gravitational waves have been derived in Ref. ||19|| . These can be obtained 
from ([39|) by specifying Ag = a, fig = k and 0g = — 7r/4. 

An alternative derivation of these equations is given in Appendix A, where 
we use the fact that the mode functions in the Heisenberg picture can be 
expressed in terms of the Schrodinger picture variables i?g, <I>g and 0g. We 
then show that the Hamilton equations for the mode functions reduce to 
those in Eq. ©. 

5 Application to a Simple Inflationary model 

Our aim is to study the growth of cosmological perturbations in the squeeze 
state formalism for a simple inflationary model. This section is mostly con- 
cerned with studying the solutions to Eq. (^). In general, when Ag, fig 
and y?g are some complicated functions of time, it is not possible to solve 
Eq. (|39|) analytically. However, before we proceed to a discussion of the nu- 
merical solution, we can get some insight into the dynamics of the system 
using analytical techniques. 

5.1 Analytic Approach 

We assume that Ag, fig and <£>g are slowly varying functions of time, i. e. for 
kr] < 1 we have Afig/fig, AAg/Ag, A<^g/Vg < 1. 

In the strong coupling or squeeze dominated regime (Ag > fig), the squeeze 
angle $g and the rotation angle 0g approach a stable fixed point (freeze 
out). The squeeze factor grows monotonically with time, which reflects the 
fact that in the course of evolution the growing mode becomes more and 
more dominant over the decaying mode. In the weak coupling regime (Ag < 
fig), the solution is oscillatory, with the squeeze factor remaining essentially 
constant and the dominant features are the oscillations of the squeezed state, 





-fig + -^(tanh % + coth%) sin2(^g - $g) 
fig - Agtanhi?gsin2(vjg - $g) . 



(39) 
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which are revealed physically as the pressure oscillations in the hydro dynamic 
fluid. 

(a) Strong coupling : Freeze out 

For A^ > there is a fixed point (in $^ and QA to the equations of 
motion 

$i' = 0* = , RL = A 5 cos 2{ n - 9$ (40) 

with 

sm2(<p e - $|) = ^(tanh^ + coth^)- 1 frign(J2 s )(£) . (41) 

^ k k 

Using this condition, we can now integrate Eq. (fiD]) for the squeeze factor to 
obtain 

J2 5 ~ /(A|-0|)V>, (42) 

so that i?g grows monotonically. Most of the squeezing occurs in the strong 
coupling regime. 

(b) Weak coupling : Oscillations 

For A^ <C Q% and taking A^, and y?g constant, we get the solution 

tan($£ - y?g) = cosa^tanf— ^(77 - r] ) + aig] - tana^ 

in which sina^ = A^/fig. In the case where a^. <C 1, this solution reduces to 
the form 

Rk = 4+^ sin2 %^-^) = 4+^|sin2(^-^). (43) 

We can consider this oscillatory solution as a reasonable approximation for 
modes well within the horizon in both the inflationary era, when A^/fig ~ 
l/k\ri\ <C 1, and the radiation dominated era, when A^/fi^ ~ 1/2. For 
these modes is constant on average, i.e. there is no net squeezing and 
perturbations do not grow. 
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For modes that cross the horizon during the matter era, where fi^ ~ Xt, 
we cannot apply this simple analysis. 

(c) An exact solution: The Bunch-Davies vacuum 

In the exponentially expanding de Sitter stage, when = k, A^ = l/\r)\ 



and ifj: = —n/2 there is an exact solution to the equations of motion (|39|) 



Rr = sinh 1 



2krj 



71 1 

<3>r = arctan ■ 



4 2krj 

Or = kn + arctan — — . (44) 
k 2kr) 



This solution corresponds to the Bunch-Davies vacuum p0| , which is an 
attractor. If the initial state (for the modes within the horizon) is not already 
highly squeezed, one finds that as the modes get driven to superhorizon scales 
they evolve toward the Bunch-Davies vacuum. In the language of squeezed 
state parameters this corresponds to the freeze-out of <3>£ and G^; we see this 
behavior in the limit k\r)\ <C 1 of Eq. 



5.2 Squeezing in a Simple Inflationary Model 

Having established that most of growth occurs on superhorizon scales, we now 
use a simple model to estimate the amount of squeezing in the perturbation 
field. We have found that all of the relevant squeezing occurs on superhorizon 
scales, i.e. when k\r]\ < 1. In Fig. 2 this corresponds to the interval [i]i x ,i]2 X }, 
where r\\ x = —1/k and r]2 X — 1/k (2/k) in the radiation (matter) era. 

The relevant squeezing occurs for the couplings for which Ag ^> £1% when 
A^ ~ z'/z. We can then integrate Eq. ( [42]) to obtain 

R^= d\nz. (45) 

J Vim 

During the inflationary era (superscript i), z can be approximated by z(rj) ~ 
(2/3) l / 2 a/l P (where l P = (87TG/3) 1 / 2 is the Planck length). The amount 
of squeezing is given by Ai? 1 ^ ~ ln(a^/a lx .) ~ ln(l/|^|fc). In the hy- 
drodynamical era, if a mode crosses the horizon during the radiation era, 
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then z(rj) ~ 2 1 / 2 a and AR™ d ~ \n(a2x/a>it) — l n (l/|r7i?|/c). If the mode 
crosses the horizon in the matter era for which z(rf) ~ a 3 ^ 2 2 1 / 2 we have 
ARf d ~ ln(W^)> and A ^F a " - (3/2) ln(a X2 /a cq ) ~ 3 ln(l/% q ). This 
last term is a poor approximation in the matter era; in fact ~ fig and the 
squeezing angle $g is not completely frozen resulting in a slower growth of 
In the case of gravitational waves (|45|) gives Ai?g ~ ln^x/aiz), which 
is in agreement with the result first obtained by Grishchuk and Sidorov in 

5.3 Numerical Analysis 

We shall now study numerically the evolution of perturbations in our simple 
model. It is important to point out that during the hydrodynamic era we 
are looking only at the collective field of baryonic matter and radiation (we 
are ignoring cold or dark matter, or any other field which cannot be accu- 
rately described by a single collective scalar field). In addition we ignore 
decoupling of matter and radiation. We do not expect to get results which 
agree completely with the highly refined calculations which already exist in 
the literature [EI]. However we do expect approximate agreement if we look 
solely at the baryonic and radiation sector of these simulations; in particular 
in the radiation era and on superhorizon scales. 

The evolution is given by the recursion relations (|3^) and we shall assume 
the following time dependence for the scale factor 

ai= ~ W ' (~°° < ">? < ~ ">?/?) inflationary era , 

a = | (^r) 2 + O^f) ' (VR < V < 00 ) hydrodynamical era, (46) 

where H = 7i/a = a' /a 2 is the Hubble constant during inflation, 9 and 77* 
are chosen such that a(r)R) = flj(?7fi) and a'(r//j) = a'^rjn), where conformal 
time t]r denotes the end of inflation (we assume instantaneous reheat). 

We normalize a such that a eq = 1 and we set 77* = 1 . As in |10| we assume 
z = aip' /1-C in the inflationary era. During the early radiation era, when most 
of the matter particles are relativistic, we take cj. = 1/3. We assume that 
there is a time, rj = r] re i, when matter particles become nonrelativistic. For 
77 > 77 rel we have c 2 s = (Sp /Se )s = 1/3(1 + 3a/4). (We checked that the 
choice of rj Te \ does not influence squeezing of the state.) The wave function 
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is continuous at 77 = rjn, which means that the functions R^, <3>£ and 0^ 
are continuous at tjr. The overall amplitude of the perturbations in the 
hydrodynamical era will be dependent on the amount of squeezing in the 
inflationary era, which in turn depends on reheat temperature specified by 
Vr- 

Evolution of the squeeze parameters 

Fig. 3 is a plot of the evolution of the squeeze factor as a function of 
the scale factor a. Most of the growth in the squeeze factor occurs on the 
superhorizon scales between the marks Ix and 2x. When kr] <C 1, the ana- 
lytic result discussed subsequent to Eq. fl4~5[) is an excellent approximation: 
i?£ ~ hxa/a\ x . 

What about subhorizon scales (k\r)\ > 1)? In inflation (when A^ <C Clt) 
R^ ~ 0, while in the radiation era R^ oscillates (see Eq. (H3|)). For the modes 
which enter the horizon in the matter era (case fcr/ eq = 0.1 on Fig. 2), the 
squeeze factor R? continues growing as Ai?g ~ Cglna, where ~ 1.3 for 
kfjeq = 0.1 and Ct — > 1.5 for kr] ccl <C 1. This means that, as a consequence of 
large coupling in the matter era (Qr < Ar), the squeeze angle remains frozen 
($r ~ const.) and the squeezing continues. Physically, this is related to the 
classical process of gravitational collapse. 

There is a critical wave vector, fc cr it?7 — 2 (which corresponds to the scale 
Aphys — 7rA eq ~ 40 (fi /i 2 )~ 1 Mpc; A eq ~ 13 (fi /i 2 )~ 1 Mpc is scale today 
corresponding to the horizon size at equal matter and radiation, Q is the 
fraction of the critical density today and h is the present Hubble parameter 
in units of lOOkm/s/Mpc). For k > k CT i t the state oscillates and R^ doesn't 
grow, and for k < k CT[t the state is frozen and R^ grows. 

Fig. 4 shows the evolution of the squeeze factor with respect to the 
squeeze phase. On subhorizon scales in inflation <3>£ grows, while R^ ~ 0. 
At the horizon crossing, the squeeze angle freezes out: $*• = ni\ (n G Z) 
becomes an attractor and, as we can read of from the figure: $i = (— 2n) 
if k — 1 (10) which is in agreement with Eq. (flop. For the subcritical case 
krj eq = 1, after the mode crosses the horizon at 2x, the angle $ ~ $* remains 
frozen and R^ continues growing. On the other hand, for kr] eq = 10 > k CT i t r} eq , 
after time 772a; the mode starts oscillating with Ai?g ~ 1. The amplitude of 
oscillations in R^ is slightly bigger than predicted by the simple formula (|43| ) ; 
the reason being that the condition for validity of Eq. (|43|) (X^/fl^ <C 1) is 
not strictly satisfied (here we have A^/f^ — > 1/2). During one oscillation 
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A$£ = 7r and remains constant on average. Looking back at Eqs. (|30|) 
and (|43|) we observe that if grows (oscillations), the squeeze operator 
produces and destroys, on average, equal number of particle pairs, i. e. there 
is no net squeezing. 

Evolution of physical quantities 

We are interested in looking at physical quantities in the hydrodynamical 
era, typically the Bardeen variable <J> B (which corresponds to the Newtonian 
potential inside the horizon) and the energy density perturbations 5e/e. In 
the standard notation 

7 = / 7^ V* ( 47 > 

e J {2tt ) 2 

and 

|<y 2 = (48) 
r> /3, H 2 -H'l 

where Zp = (87rG/3)~5 is the Planck length. 

To make contact with the existing work on power spectra from inflation, 
we plot in Fig. 5 the growth of the power spectrum \8^\ 2 defined in Eqs. ( (4"7|) 
and ( fi8|) against the scale factor for the modes: krj eq = 0.1 and kr] C(l = 3. 
On superhorizon scales, during the radiation era (In a < 0), the power grows 
as 1 6%\ 2 oc a 4 oc r/ 4 , which agrees with the estimates based on Eq. (|45|). In 
the matter era, for the modes kr] eq < 2, the power grows as \5j:\ 2 oc a 2 oc r/ 4 , 
while for kr] eci > 2, the state start oscillating and the growth becomes very 
slow. 

Fig. 6 shows the power spectrum at two different time slices: rj = O.lr^q 
and rj = 0.5r] eq . The spectrum is scale invariant, |5^| 2 ~ k, on superhorizon 
scales and \5^\ 2 ~ k~ x on sub horizon scales. The turning point, caused by 
the oscillations of the squeezed state (see Fig. 5), is at ki] eq ~ 4 for both 
time slices. The first dip in the power spectrum is at krj eq ~ 8 — 9, which 
corresponds to the wavelength A ~ (0.8 — 0.9)A eq — 11 (fio^ 2 ) _1 Mpc. These 
dips correspond to the acoustic oscillations in the fluid. 
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6 Comparison with previous work 



The features of the power spectrum just discussed are those expected. We 
obtained the correct growth on superhorizon scales and found acoustic oscil- 
lations in the modes which reenter the horizon in the radiation dominated 



era, as described, for example, Bardeen et al. in Ref. [l[ and in Ref. pl| . We 
have simply illustrated that these phenomena can be described in a different 
way in the squeezed state framework. 

For a more direct comparison with the work of Grishchuk and Sidorov in 
Ref. @ , in particular their discussion of "desqueezing" , we treat analytically 
a model in which matter and radiation instantaneously decouple. We work 
with the action of Eq. ( p0|) as in 0. We take 

a = V < V2 < (inflation) (50) 

j ~~V2 < f] < r/i (radiation) (51) 



rj 2 



" (jjtji)*' V> 2 Vi (matter). (52) 

The most convenient way to solve explicitly for the squeeze factor R is to solve 
the equation for v derived from the action and then to use the transformations 
relating the two sets of variables derived in Appendix A. The solutions for v 
in the three eras are 



V; 



U l - i>~ ikr> (53) 



v r = _ V / S e ~^ (54) 

and 7r = v'. The normalizations are chosen so that v'v* — v'*v = —% in each 
case. Matching v and v' ( and therefore R and $) continuously at each of the 
two boundaries we obtain the following expressions for the squeeze factor to 
leading order in kr\2- 

sinh2 ^= ' ^<V2 (56) 

1 = 4(fcr/ 2 ) 4 V 2 ~ C ° S 73 J ' ~ V2 < V < Vl ^ 



sinh 2 Rr = — - -r l (a + P)( (— ) 4 + 7r^(— 



2 
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, ( y 

+ 7 W 



) + (a - 8 K0 2 + i(^ ( T )1 )}'' >2 ' fe(58) 



where 



a 



1 

12 



^ 5 + 



7 



1 + 



2(kr ll 



\{kr]i 



+ 1 + 



2(% 



cos 



COS 



COS 



2% 

7T 



4fc?7i 2fcr/i 
— pr- sin 

- — = sin — =- 

\/3 V3 



sin 



v/3 
2fcr/i 



(59) 
(60) 
(61) 



From the second expression in ([58) we see that the squeeze factor is 
modulated only slightly in the radiation era in agreement with what we 
found earlier. For the matter era however one can show that the coefficients 
a + (3 and 7 of the terms which grow with 77 vanish when the condition 



km 2km 
— = + arctan — ■=• 



in 



integer) 



(62) 



is satisfied. This leads to a significant amount of "desqueezing" of these 
modes as the squeeze factor for these modes is given approximately by 



sinh 2 R t = sinh 2 R° r { ^ 

V 



2 1 + 



(kri) 2 



1 + 



(63) 



(2fcT?l) 2 



where R~ = R^{2rji) is the squeeze factor at the decoupling. In terms of the 
scale factor, 



R° t - -ln( 
k 2 



for kr] ^> 1. 



(64) 



The existence of this "desqueezing" is again a familiar phenomenon ex- 
pressed in a different set of words. When one matches the oscillating solutions 
of the radiation era onto the growing and decaying solutions of the matter era 
one finds that certain modes match completely onto the decaying solution. 
In fact this is the simplest way to derive the condition (|62|) above. These 
modes lose power and we have approximate zeros in the power spectrum. 
These oscillations in the power spectrum are known as Sakharov oscillations 
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TQ. In order to obtain the position of the zeroes, we solve Eq. (B^) and 
obtain k7] iec = 2krj! = {6.36, 16.7, 27.4, 38.2, 49.1, 59.9, 70.8, ..}, which cor- 
respond to today's scales: A = {89, 34, 20.7, 14.8, 11.5, 9.5, 8.0, ..} h^Mpc 
(h is the Hubble constant today in units lOOkm/s/Mpc ). The occurrence of 
these oscillations depends crucially on the matching at the inflation-radiation 
transition. In order to match purely onto the decaying solution (in the mat- 
ter era) , one must have standing wave solutions in the radiation era and this 
in turn depends on having the correct input from the inflationary epoch. It 
is indeed the squeezing of all of the physical momentum out of the super- 
horizon modes during inflation that produces the standing waves at the end 
of inflation, which one requires to produce this effect. 

Grishchuk and Sidorov suppose this crucial ingredient to be missing in 
standard treatments of the growth of perturbations. They claim that in- 
correct assumptions about the perturbations produced by inflation are often 
made which lead to traveling wave solutions in the radiation dominated era 
and the resultant absence of these Sakharov oscillations in the power spec- 
trum. For example, Grishchuk states in |6| that "the unavoidable property 
of squeezing manifests itself in the fact that the phases of primordial den- 
sity perturbations are fixed and correlated, in contrast to the usually made 
assumption that the phases are distributed randomly and evenly. In other 
words, the primordial density perturbations, similarly to the case of grav- 
itational waves, must form a set of standing waves with definite phases". 
[Our italics]. In fact these two points are not in conflict. Indeed there are 
"standing waves with definite phases" , but there are other phases which are 
distributed randomly and evenly. One must be careful about which phases 
one is talking about. Each standing wave has a "phase of oscillation" which 
distinguishes among solutions which are at different points in their period of 
oscillation. This is the phase which is fixed (relative to the time of horizon 
crossing) in inflationary cosmologies. 

However, inflation does not predict the location of the nodes in the stand- 
ing wave. There is another "spatial" phase which distinguishes among stand- 
ing waves which differ by a translation in space. Since the wavefunction as- 
signs equal probability to solutions which differ only by a translation, one 
can choose a random spatial phase. This amounts to making a particular 
random choice of S(x) from among the many possible ones. 

We are aware of one place in the literature where an error is made regard- 
ing which phases are random. In a passage in (preceding the paragraph 
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containing Eq. (7)) Peebles argues that the temporal phase of the standing 
waves may be taken to be random. This statement is incorrect and, to the 
extent that Grishchuk's criticisms refer to it, we are in agreement with him. ^ 
However, this is an isolated error and is not of significance in either the work 
of this author or others who produce detailed predictions based on specific 



models (see e.g. Ref. [21 



Typically the correct standing wave solutions are used without making 
reference to the squeezed state terminology. That this so can most simply 
be seen by the fact that the usual Bunch-Davies vacuum matched onto the 
oscillating radiation era solutions gives precisely the standing waves noted 
by Grishchuk and Sidorov. In |[21|| , for example, the matching is described in 
terms of growing and decaying modes in the radiation era, but amounts to 
the choice of standing waves and indeed both the acoustic oscillations and 
Sakharov oscillations which result are seen in these simulations. The reason 
why so little attention is paid to these features is that they occur only in the 
baryonic component of matter and are almost completely swamped in dark 
matter dominated models. It is an interesting possibility that this difference 
might be exploited to distinguish between baryonic and dark matter domi- 
nated models. Attempts have in fact been made to look for these Sakharov 



oscillations but the results are inconclusive 24 



The other important claim of Grishchuk and Sidorov is that these fea- 
tures can be said to be of a distinctly quantum mechanical origin. Speaking 
of desqueezing, they state in @ that "we relate this quantum effect to the 
effect of the so-called Sakharov oscillations known in the classical theory of 
matter-density perturbations" . In || Grishchuk opines that "it is quite pos- 
sible that the very specific properties of the large scale density perturbations 
related to their quantum mechanical origin can be revealed in the appropriate 
observations" . [Our italics] . We will attempt to clarify this question in the 
next section. 



2 We are grateful to Jim Peebles for a discussion of this point. 
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7 The Classicality of Squeezed States 



A squeezed state seems to be an especially quantum mechanical state. It is 
not well localized in p and q and therefore cannot be represented by a point 
in classical phase space. It may instead be viewed as a coherent superposition 
of many localized wave packets. It is very unlike the archetype classical state 
- the coherent state — being very squeezed in one variable. It is this feature 
which generates so much interest in these states in quantum optics and other 
areas of physics and leads to their characterization as very "non-classical" 

PH. 



7.1 Quantum coherence 

An important question is: How can the quantum coherence of the squeezed 
state manifest itself in physical processes? To clarify this question let us make 
a few general remarks about quantum coherence. A wavefunction ip(x,t) 
assigns probability z^* (a:, t)^(x, t) to the states \x). If all one ever asked 
about were the probabilities assigned to states \x) at time t, one would be 
working with the equivalent of a classical probability distribution. Quantum 
coherence comes into play when one asks, for example, about probabilities 
assigned at time t to states other than x eigenstates. At this point know- 
ing tp*(x,t)ip(x,t) (the probability distribution in x space) is not enough, 
and one needs the information provided by the complex function ip(x,t) to 
generate new probabilities. One can say that this is because the state is a 
"coherent superposition" of x eigenstates. It is possible to put the system 
in an incoherent superposition of x eigenstates by representing its state as a 
density matrix of the form \x)p(x)(x\. In this case the probabilities p(x) are 
all you need to know. 

In the case of the coherent superposition if)(x), one can avoid all ques- 
tion of quantum coherence by limiting one's attention to the probabilities 
assigned to the x eigenstates at time t. However, the nature of the time 
evolution can make the quantum coherence hard to avoid. A well known 
example is the double slit experiment. If one starts at some time t\ knowing 
ip(x,ti) for the electron before it passes through the slits, then at a late 
time t 2 , ip*(x,ti)ip(x,ti) will always give the probability assigned to what- 
ever the initial state \x, ti) evolves into under time evolution to £2. One 
can avoid questions of quantum coherence at time t 2 by limiting one's at- 
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tention to these evolved states. The problem is that in the double slit ex- 
periment the \x) states evolve into something very complicated, and one's 
attention (eg measurement) is focused on other simpler states (such as states 
which are eigenstates of x at t 2 ). One thus requires a knowledge not just of 
if}*(x,ti)ip(x,ti), but of the full complex phase information in ip(x,ti). This 
is when quantum coherence is important. 

However, when the evolution of each basis state (with respect to which 
the initial wavefunction is expanded) is simple enough, one can realistically 
expect to limit one's attention to whatever states this initial basis evolves 
into. This allows one to regard the square of the wavefunction as giving a 
classical probability distribution, and avoid any question of quantum coher- 
ence. A particular example of this simple evolution is when the state is WKB 
classical. 



7.2 WKB classicality of squeezed states 

Consider the q representation of the squeezed state in the static inverted 
harmonic oscillator which we considered earlier 

^j(q) = Ne~^ B+iC)q i, (65) 

where 

N = (^-) 4 , B — — ) — , C = tanh 2r. (66) 

\hnj ' cosh2r' v ; 

We will show that for large squeezing this wave function is very classical 
in the WKB sense and becomes increasingly so with time. The wavefunction 
can be written 

if,(q) = p (q)e iS ^. (67) 

If S(q) varies much more rapidly with q than p(q) the state is a WKB state 
for which 

p\^} ~ (hd q S(q)m. (68) 

To the extent that this holds the state assigns momentum and position si- 
multaneously according to 

p(q)=hd q S(q). (69) 
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While p(q) need not be localized, it does represent a distribution in classical 
phase space which evolves classically in the WKB limit. 
For the evolved jstate given by Eq. (pSI) we have 



p(q) = Ne~ B ™ (70) 

S(9) = ("I 

The WKB condition is met when the quantity p(d q S(q)/d q p(q)) is large. 
From Eq. flB5p we find 



d q S(q) 



d q p(q) 



C 

— = sinh2r. (72) 
B 



Therefore as the initial state evolves and becomes more squeezed, it also 
becomes more classical in the WKB sense. 

Equivalently this can be seen from Eqs. (§) - ( |TT| ) since they imply 

\{V{t)\qp\V{t))\ = |(l + sinh 2 2r)5 ~ ^ e 2t . (73) 

This just expresses more directly the effective irrelevance of the noncom- 
mutativity of the position and momentum operators on the state for large 
squeezing. 

It is precisely these properties of the inverted harmonic oscillator which 



were used by Guth and Pi in [ I !j to illustrate how a quantum mechanical 
state can be treated in certain cases as an ensemble of classical states. This 
WKB classicality means that the squeezed state can be approximated in its 
evolution as a classical phase space distribution, as long as one only measures 
classical quantities. When a particle in a spread out WKB state interacts 
with another system which responds to (or "measures" ) the value of p or q, 
one can predict the outcome using only the probability distribution in classi- 
cal phase space. One does not need to know the complex phase information 
contained in the full wavefunction, so questions of quantum coherence do not 
arise. 

In fact, it is well established that when such a measurement takes place 
correlations are set up which cause the quantum coherence to be lost (see for 
example PS[ ) . From that point on the particle is in a density matrix rather 
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than a pure state, and the possibility of observing the effects of quantum 
coherence is even more remote. 

In quantum optics, where squeezed state of the electromagnetic field can 
be produced, one can not think in terms of a classical probability distribution. 
This is because the electromagnetic fields are measured by the absorption of 
photons by atomic systems. These interactions typically do not amount to a 
measurement of the classical field variables, and so quantum coherence effects 
are observed. 

The crucial question then is: when matter interacts with a density field 
in a squeezed quantum state, does it respond to (or measure) the classical 
field values or something else? We have given this question some thought, 
and find it hard to see anything other than very classical processes in these 
interactions. The matter, after all, evolves according to the values of things 
such as the Newtonian potential, which is local in the field variable. Fur- 
thermore as the universe evolves the matter responds to the perturbations, 
correlations will be set up which destroy the initial coherence as discussed 
above. 

If one wishes to show that the initial quantum coherence of the squeezed 
state is of physical importance, one must demonstrate interactions which 
measure something other than classical quantities before the ordinary inter- 
actions destroy the quantum coherence. It would be very interesting if this 
could be done, but we do not see how. 

The particular features of the power spectrum discussed by Grishchuk and 
Sidorov are not the result of quantum coherence. They are features which 
appear in individual classical solutions (eg properties of each trajectory in 
classical phase space) and do not represent quantum interference among dif- 
ferent classical solutions. The physical origin of the fluctuations (the vacuum 
fluctuations) is quantum mechanical but their known physical effects are in- 
distinguishable from fluctuations from a classical stochastic field. 

Regarding the phases of modes which oscillate inside the horizon, these 
are predicted regardless of whether there is quantum coherence. The predic- 
tion is based on the fact that the modes in question have spent a long time 
outside the horizon, where there is one growing and one decaying solution 
to the equations of motion. The growing component becomes completely 
dominant for the modes which are amplified during inflation. This growing 
solution has a uniquely determined oscillatory behavior when it enters the 
horizon, and thus the phase of the oscillations is predicted. The original work 
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on this subject has correctly accounted for these predictions PJ. 

The quantum squeezing is also a consequence of the presence of one grow- 
ing and one decaying solution, but that does not mean that observing the 
phases of the oscillatory behavior amounts to a test of quantum coherence. 
An incoherent superposition (such as would result from the establishment 
of correlations with particles and photons mentioned above) would provide 
the same results, as long as each mode was dominated by the growing solu- 
tion. The particular features of the power spectrum discussed by Grishchuk 
and Sidorov are only quantum mechanical in origin in the mundane sense 
in which all perturbations in inflation are. The physical origin of the fluc- 
tuations is quantum mechanical but they are in their known physical effects 
indistinguishable from fluctuations from a classical stochastic field. 

8 Conclusion 

We developed the squeeze state formalism to study the growth of cosmolog- 
ical perturbations. The formalism is then applied to a simple inflationary 
model with baryonic matter. We discussed how the standard features, such 
as acoustic oscillations and Sakharov oscillations, are characterized in the 
squeeze state formalism. At late times density perturbations are semiclassi- 
cal and — for all practical purposes — can be well represented by a classical 
probability distribution function. 

Confusion can be avoided if one keeps in mind that there are three very 
different phases which enter into the discussion. Firstly, there is the complex 
phase of the wavefunction. To the extent that the system being studied 
behaves like a classical probability distribution, this phase can be ignored. 
Secondly, there is the phase of oscillation of standing waves in the density 
field. These are very precisely fixed in inflationary cosmologies and this 
can lead to predictable Sakharov oscillations at late times. Note that this 
second phase is a classical phase. Thirdly, there are classical phases for each 
Fourier mode which correspond to translations in physical space. Since the 
inflationary universe assigns equal probability to density fields which differ 
only by a translation, these spatial phases are random within the linear 
approximation. 

The use of squeezed states in a cosmological setting was first advocated 
and implemented by Grishchuk and Sidorov to calculate the power spectrum 
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of primordial gravitational waves [|T^]. The treatment is entirely analogous 
to that of cosmological perturbations; it is possible to reduce the problem 
to, again, quantizing a scalar field with a time dependent mass (z = a). 
The power spectrum of this scalar field exhibits oscillations on certain scales. 
It is possible, as Grishchuk claims, to predict the position of the dips in 
the power spectrum. However this feature is also present in the standard 



Heisenberg formalism as treated by Abbott and Harari [26]. The squeezed 



state formalism gives us an intuitive way of looking at the generation and 
evolution of cosmological perturbations. However the formalism we have 
developed is not restricted to cosmological applications; the equations of 
motion that we have derived are quite generic of systems with quadratic 
hamiltonians that can be put in the form of Eq. . 
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A Relating the Heisenberg and Schrodinger 
pictures 

In this appendix we show how to parametrize the Schrodinger picture vari- 
ables in terms of the squeeze state parameters R^, and 0r. We then 
demonstrate that the classical equations of motion for the mode functions re- 
duce to the evolution equations for R^, and 0g derived in Sec. 3 (Eq. (|39l)). 
This shows how the Schrodinger picture problem can be reduced to solving 
the classical equations of motion with an appropriate reparametrization. 
The Heisenberg picture operators v{x,rj) and 7c(x,r]) can be written as 

v{x, rj) = W(r), r} )v{x, r] )U(ri, rj ) = J -^m^ {ut{r,)a^ + U-fcfa) !^) 
7c{x,r]) =Wt(77, 77o)7r(x, r)o)U(r), Vo) = f + w _ % {ri)(2_£U) 

It is now easy to show, using the Heisenberg equations of motion for v and tt, 
that the mode functions uAf]), wAjf) satisfy the following Hamilton equations 

z' 

w k + ~ u k 

- C l k2u k- ~ W k- ( 75 ) 

These are the configuration and momentum variables of the classical field 
theory given by the action in Eq. ([H]). With the initial choice of uA'Hq) = 
(2/c) -1 / 2 and wAt]q) = i(k/2) 1 ^ 2 , corresponding to an initial (right moving) 
traveling wave, the solution to Eq. (|75[) are uniquely defined for all times. At 
V = f]o we obtain the Schrodinger picture operators At some later time 
r) we have 

=-i>/fk('?)-als(^)]. ( 76 ) 

where a^rj) and a) Jrf) are the Heisenberg picture annihilation and creation 
operators defined by 

a R (rj)=U\rj,rj ) H U(rj,rjo) = 7^(%)<St(% %K<S(% ^)W(e s ) 

= cosh R z e- i@ *a % - sinh R n e^ +2 ^ a^77) 
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From Eq. ([74]) we then get 



Comparing Eq. ( ff"4l 



g (cosh %e" ie E - sinh R % e~ i( - e ^ +2 ^ 
+al f . (cosh i^e ie s - sinh R^ e i+ 2 ^ 



HI 2 



-k 



% % (cosh%e" ie s + sinh%e-^ e s +2 *s) 
(cosh R R e ie K + sinh R^®^ 2 ^)] . (78) 
with Eq. fl78D we can identify the mode functions to be 
= -j=(coshR n e ie k - sinh flgeW" 2 **) 



( cosh Rre ie k + sinh R t e i{ - @ % +2 *%- 



(79) 
(80) 



and these define the transformation that we seek between the Schrodinger 
picture variables and the Heisenberg picture mode functions. It is now a 
matter of algebra to show that Hamilton's equations for the mode functions 
( |75|) give the equations of motion for R^, $g and 0£ (|39|). 



B Invariance of the equations of motion for 

R % , * £ and 9 £ 

Here we show that for the two actions fl20|) and (|2T|) differing by the total 
derivative term ((z'/z)v 2 )', the equations of motion for R^, $^ and have 



invariant form. For the action (|2T| ) A^, f2g and ip^ are defined in Eq. (|2 
On the other hand for the action (^) we have 



2kz 



71 

~4 ' 



(81) 



Even though canonically related Hamiltonians give different evolution for R^, 
and G^, the physically measurable quantities are invariant. We have not 
investigated how generic is the invariance of the equations of motion (|39|). 
We leave this as an exercise to an inquisitive reader. 
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Figure captions 



Fig. 1. 

Phase space trajectories for a classical upside-down harmonic oscillator. 
The presence of one growing and one decaying solution produces a "squeez- 
ing" effect even at the classical level. The circular region shown evolves with 
time into the squeezed shape above it. 

Fig. 2. 

Evolution of scales in an inflationary universe model. x c denotes the 
comoving scale and a# the end of the inflationary stage. The perturbation 
is on super horizon scales in the interval [a(rj x i,r] x2 ]. 

Fig. 3. 

Evolution of the squeeze factor R as a function of the scale factor log a in 
an inflationary universe model for two scales: krj eq = 0.1 and ki] eq = 3. Most 
of the growth occurs on superhorizon scales (period between the marks lx 
and 2x). 

Fig. 4. 

The squeeze phase vs. squeeze factor ($ — R) diagram for two scales: 
kr] eq = 1 and kt] eq = 10. The squeeze angle freezes out on superhorizon 
scales. On subhorizon scales it exhibits two types of behavior: for scales 
bellow critical (k cr i t r] eq ~ 2), the $ — R curve exhibits oscillatory behavior, 
and for scales above critical the phase remains frozen. 

Fig. 5. 

The growth of the power spectrum \5^\ 2 against the scale factor a for 
/c?7 eq = 0.1 and kr] eq = 3. In both cases we observe the same power law 
growth on superhorizon scales. In the subcritical case (ki] cq = 0.1) the growth 
continues (with somewhat slower rate), while in the supercritical case (krj eq = 
3) the power exhibits oscillations after the horizon crossing. 

Fig. 6. 

The snap-shot of the power spectrum (log \5^\ 2 - log k plot) for two times: 
a = 0.1a eq and a = 0.6a eq . The spectrum is scale invariant on superhorizon 
scales: \5^\ 2 ~ k and after the turning point at kr/ cq ~ 4 it shows oscillatory 
behavior and decays as \5^\ 2 ~ k' 1 . 
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